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Abstract 

We demonstrate the correspondence between quantum mechanical and quantum 
field theoretical descriptions of Mossbauer neutrino oscillations. First, we compute 
the combined rate V of Mossbauer neutrino emission, propagation, and detection in 
quantum field theory, treating the neutrino as an internal line of a tree level Feynman 
diagram. We include explicitly the effect of homogeneous line broadening due to 
fluctuating electromagnetic fields in the source and detector crystals and show that 
the resulting formula for F is identical to the one obtained previously [l| for the case of 
inhomogeneous line broadening. We then proceed to a quantum mechanical treatment 
of Mossbauer neutrinos and show that the oscillation, coherence, and resonance terms 
from the field theoretical result can be reproduced if the neutrino is described as a 
superposition of Lorentz-shaped wave packet with appropriately chosen energies and 
widths. On the other hand, the emission rate and the detection cross section, including 
localization and Lamb-Mossbauer terms, cannot be predicted in quantum mechanics 
and have to be put in by hand. 



1 Introduction 



The possibiUty of exploiting the Mossbauer effect in weak interactions to enhance the small 
neutrino cross sections S [si S, 0, [gI has recently received considerable interest, both from 
the experimental side p, la, l7|, l8|] and from the theoretical side Q, i, 0, EH, Q [lal, E3, H [lej • 
In the proposed experiment, neutrinos are emitted from atoms embedded into a metal 
crystal and absorbed by ^He atoms embedded into a similar crystal. With very optimistic as- 
sumptions on the source activity (1 MCi), the fraction of recoilfree emissions and absorptions 
(0.28 each), and the achievable spectral line width (AE/E ~ 10^" eV/18.6 keV ~ 5-10-^^), 
it has been estimates that an event rate of 10'^ per daycould be achieved for a detector 
containing 1 g of ^He and placed at a baseline L = 10 m pill These events could be counted 
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^It has been suggested recently that it might even be possible to reach a line width of 0(10"^^ eV), 
corresponding to the natural line width of tritium decay 17, 18, [l9| . This would imply an additional 
enhancement of the event rate by a factor of lO^f, allowing for smaller sources and detectors, or for longer 
baselines. However, the arguments given in 17, 3, [3l in favor of this additional enhancement have been 
disproven in ref. [20]. 



by observing the subsequent decays of the produced '^H in the detector, or by chemically 
extracting and counting the number of produced atoms. However, it is far from clear 
whether the above experimental performance can be achieved in practice, and, in fact, the 
event rate may well be many orders of magnitude smaller [sl, so that the question of whether 
a Mossbauer neutrino experiment can be realized in practice is still open. 

In spite of this, Mossbauer neutrinos have already now proven to be an excellent test 
case for studying the quantum mechanics and quantum field theory of neutrino oscillations 
theoretically. In particular, their extremely small energy spread of O{10-^^ eV) 0, [2lj has 
led to the question whether a coherent emission and absorption of different neutrino mass 



eigenstates, which is a prerequisite for oscillations, is possible [ll|, [IJ]. Even though a de- 
tailed quantum field theoretical treatment, requiring no a priori assumptions on the neutrino 
wave function, shows that oscillations do occur in a Mossbauer neutrino experiment (ll. [l3|. 
such an approach also reveals that Mossbauer neutrinos are special because many of the 
assumptions and approximations that are commonly made in the theoretical treatment of 
conventional neutrino oscillation experiments are invalid for them. 

In this paper, we will use the example of Mossbauer neutrinos to discuss the correspon- 
dence between quantum mechanical and quantum field theoretical approaches to neutrino 
oscillations. In sec. [2], we will derive the combined rate of Mossbauer neutrino emission, 
propagation, and detection in quantum field theory (QFT). We will for the first time explic- 
itly include homogeneous line broadening effects arising from fluctuating electromagnetic 
fields in the solid state crystals forming the Mossbauer source and detector, and we will 
show that, as anticipated in ref. [H, the result is identical to the one obtained for inhomoge- 
neous line broadening due to crystal imperfections. We will then derive the same result from 
quantum mechanics (QM) in sec. [31 treating the neutrino as a wave packet. A comparison 
of the QFT and QM approaches will show that QM is inferior to QFT because more ad 
hoc assumptions are required, e.g. on the shape and width of the neutrino wave packets. If, 
however, all parameters are chosen appropriately in the QM formalism, the QFT result can 
be reproduced. In sec. HI we will discuss our results and conclude. 



2 Mossbauer neutrinos in quantum field theory and homogeneous 
fine broadening 

To compute the amplitude for Mossbauer neutrino production, propagation, and detection 
in QFT, we follow the formalism developed in [l| and consider the Feynman diagram shown 
in fig. [TJ Here, the external lines correspond to the ^H and '^He atoms in the source (S) and 
the detector (D), while the internal line describes the propagating antineutrino. Since we 
are mainly interested in the phenomenology of Mossbauer neutrino oscillations, we avoid 
an explicit treatment of solid state binding forces and instead assume the external particles 
to reside in the ground states of simple harmonic oscillator potentials, with oscillator fre- 
quencies of the order of the Debye temperature ~ 600 K ~ 0.05 eV of the respective 
crystals [sl, 0]. It is known from the theory of the classical photon Mossbauer effect 22 



that this model provides qualitatively correct results, even though it is, of course, insufficient 
for computing a precise prediction of the total event rate. If we denote the masses of the 
external particles by {A = {H, He}), their average positions by {B = {S,D}), the 
harmonic oscillator frequencies by lja,b, and the ground state energies by Ea,b, the wave 
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Figure 1: Feynman diagram for neutrino emission and absorption in the H- He system 



functions corresponding to the external legs in fig. [T] are given by, 
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Due to interactions of the atoms with their surroundings, Ea,b will not be constant in 
time, but will fluctuate around the average zero point energy Ea,b,o = + \^a,b [21 
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These fluctuations are, for example, induced by random thermal spin 
flips of neighboring atoms. They are generally referred to as homogeneous line broadening 
effects because, as we will see below, they limit the achievable sharpness of the Mossbauer 



resonance. To describe homogeneous line broadening, we make the replacement [23 
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in eq. ([T]), where 
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is integrated phase shift induced by the fluctuations of Ea^b- Note that this approach 
accounts only for homogeneous line broadening due to solid-state effects, but not for broad- 
ening due the natural line width. The latter effect (which is theoretically interesting, but 
completely negligible in the ^H-^He system) has been studied in detail in ref. The 



3 



transition amplitude corresponding to fig. [H including the modulation factors ([3]), is 
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Here, are the neutrino mass eigenvalues, Ugj are elements of the leptonic mixing matrix, 
and the nonrelativistic (i.e. momentum-independent) electron spinors are denoted by Ue^s 
for the electron that is emitted in decay in the source, and by Ue,D for the electron that 
is destroyed in the neutrino capture process in the detector. The matrix elements A4g and 
A^^ are given by 
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where Gp is the Fermi constant, 6c the Cabibbo angle, and u^ b (with A = {H, He}, 
B = {S,D} as before) are the non-relativistic ^H and ^He spinors. The vector and axial 
vector (or Fermi and Gamow- Teller) nuclear matrix elements are My = 1 and M^. ~ v^, 
respectively [23, 28|, and the axial-vector coupling constant is qa — 1.25. The quantity 
ijjeiR) gives the value of the anti-symmetrized atomic wave function of ^He at the surface 
of the nucleus, while the factor 
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accounts for the fact that the spectator electron in bound state ^H decay and induced 
orbital electron capture on ^He changes from the Is state of ^H into the Is state of "^He, or 
vice- versa. 

The spatial integrals in (jlj) yield a factor exp[— p^/2(Tp] exp[zpL], with the effective mo- 
mentum uncertainty (Xp of the experiment defined by 
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and with the basehne vector 



To evaluate the three-momentum integral, we employ the Grimus-Stockinger theorem 
which states that, for any three times continuously differentiable function '0(p) (p G M^), 



with tp and all its first and second derivatives decreasing at least as l/|pp for |p| 
following relation holds for any real number A > 0: 
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Effectively, this formula gives the form of the Feynman propagator for propagation over 
macroscopic distances. The parameter A corresponds to squared modulus of the on-shell 
momentum component of the propagating particle. Applying the Grimus-Stockinger theo- 
rem to our expression for A, we find 

—7 ^ r°° 



dpo exp 
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with pj = (po, JpI 



m?- L/L) and with the constant 
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containing the wavefunction normalization factors from eq. ([T]) and the numerical prefactors 
that have arisen in the xi and X2 integrations. Since we do not know the exact form of the 
modulation factors fA,B{t), we cannot evaluate the time integrals at this stage. However, 
ultimately, we are only interested in the transition rate F, which is proportional to (AA*), 
the statistical average of AA* over all possible and ^He states in the source and the 
detector. This expression can be simplified using statistical arguments. In particular, when 
evaluating it, we encounter the quantity 
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and a similar term from the detector-related modulation factors. Here, ti and ti are the time 
variables appearing in the expressions for A and A*, respectively. To shorten the notation, 
we have defined a quantity AEsit') = Esit')-Esfi = [^H,s(^')-^Hc,5(^')]-[^H,s,o-^He,5,o], 
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which gives the deviation of the energy of the neutrino emission hne from its mean value 
at time t' . Following 23|, we assume AEs{t') to be a Gaussian random variable centered 
around zero: 



{AEsit')) = 0. 



(13) 



Moreover, we assume fluctuations at different points in time to be uncorrelated (Markovian 
approximation), which implies 



[AEsit') AEs{t"))=jsS{t'-t"). 



(14) 



This is a good approximation if the correlation time of the fluctuations is much smaller 
than all other time scales appearing in the problem, in particular the tritium life time and 
the running time of the experiment. The constant 75 will turn out to be the width of the 
neutrino emission line. Proceeding along the lines of refs. 23|, ISOj, we expand ( fT2l) into a 
Taylor series and obtain 



Bs{h,h) = J2-^ I dt^''>---d&\AEs{t^'^)---AEs{t^-'>)) 



(15) 



One can now use the assumption that AEs{t^^^) is normally distributed around zero to show 
that the n-point correlation functions on the right hand side can, for even n, be rewritten 
by splitting them into products of two-point functions (which can be evaluated by using 
( JT^ ) and summing over all (n — l)(n — 3) ■ ■ ■ 3 ■ 1 = n!/[2"'/^(n/2)!] distinct combinations of 
such two-point functions. For odd n, the ra-point correlation functions can be transformed 
into products of [n — l)/2 two-point functions and a one-point function, which is zero by 
virtue of eq. (fT3|) . Therefore, Bs(ti,ii) takes the form 
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The analogous expression for the detector-related modulation factors is 
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the expression for (AA*) now becomes 
AT- 
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We can rewrite the squared (5-function as T/27r ■ (5(po — Po) (with T the total running time 
of the experiment), and use the remaining (5-factor to evaluate the pQ integral. We are 
left with the po integration, which receives its main contribution from the region where 
\Esfi — Po\ ^ Is and I-E'd.o — Pol ^ Id due to the Lorentzians on the right hand side of 
eq. (fT9l) . Since ^3,0 CTp and js,d ^ Esfi^Ej^ Q, the spinorial factors as well as the real 
exponential that will lead to the generalized Lamb-Mossbauer factor and to the localization 
term are almost constant over this region and may be replaced by their values at 
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7s/27r 7D/27r 



Ijk 



dp'' 



(p° - ^5,0 



exp 



' 2p0 



2iT (^Es,o - Ed,o) 



(7S+7d)^ 
4 



75 + 7D / .(S) 



^jk ) 



1 {Af - Af) [(E^,o - E^,o)(75 - 7d) ± ^ 

2 Es.-Ed.^±i^-^ 



(7S+7d)^ 



(20) 



(21) 



with the abbreviations 
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and with the oscillation and coherence lengths 
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In eq. ( I2T!) . the upper (lower) signs correspond to A.rn?-f^ > (Am|^ < 0). Thus, the 
transition rate F for a Mossbauer neutrino experiment dominated by homogeneous line 
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broadening is, according to Fermi's Golden Rule, 
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is the rate of bound state H decay, and 

Bo = AnGlcos^e, \MR)\''{\Mv\'' + 91\Ma\^) kj, 
is related to the cross section for induced orbital electron capture on free ^He by 31 

a{E^) = Bo p(-Ep,res) • 



(25) 
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Here p(-Ep,rcs) is the spectral density of incident neutrinos, i.e. the number of neutrinos 
per unit energy interval, at the resonance energy for this case, -Ep,rcs = Q + Er (where 
Q = 18.6 keV is the Q-value of the process and Eji is the recoil energy transferred to the 
atom). The quantities Ys and Yd in eq. (l24l) are given by 



Y, 



B 



for B = {S,D}. 

As anticipated, flMl) coincides precisely with the corresponding expression for the case 
of inhomogeneous line broadening, given in eq. (41) of ref. [l|jf| We find again a Breit- 
Wigner-like resonance term, which suppresses Mossbauer transitions if the central energies 
Esfi and E^ o of the emission and absorption lines differ by more than the average line 
width (7s + 7d)/2, and a factor 
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(30) 



which we interpret as a generalized Lamb-Mossbauer factor (or fraction of recoil-free emis- 
sions/absorptions), multiplied with a localization term. The latter can be neglected if a"^ ^ 



^In the present work, we have chosen to present T in a form where the Breit-Wigner term is factorized 
out of the term containing the osciUation and coherence exponentials. It is straightforward to check that 



this form is identical to the form used in ref. [ij . 
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Am|;,, or equivalently, if L°^jj. ^ AnaxE/ap, where = l/2ap is the spatial delocahzation 
of the emitting and absorbing atoms. For E = 18.6 keV and ap ~ {niuOnY^^ ~ 7 keV, it 
is clear that this inequality is easily fulfilled since is of the order of the interatomic dis- 
tance, while -Z^B j-fc ~ 20 m for oscillations driven by the atmospheric mass squared difference 
Am|^ and -^^ jfc ~ 600 m for oscillations driven by the solar mass squared difference Am|^. 

The factors A^f ^ in eq. flM|) contain the oscillation exponentials and the decoherence terms 
which describe the effect of wave packet separation due to the different group velocities asso- 
ciated with different neutrino mass eigenstates. However, it is easy to see that decoherence 
is not an issue in any realistic Mossbauer neutrino experiment because the corresponding 
coherence lengths are of O{10^^ km). 

The fact that the formula for F is identical for the cases of homogeneous and inhomoge- 
neous line broadening implies that these two situations cannot be distinguished experimen- 



tally. This confirms a more general theorem by Kiers, Nussinov, and Weiss [32] which states 
that it is impossible to distinguish an ensemble of neutrino wave packets with identical mo- 
mentum distributions from an ensemble of plane wave neutrinos whose individual momenta 
follow the same distribution. In fact, the density matrix describing the ensemble is identical 
for both cases. Applied to Mossbauer neutrinos, the case of neutrino wave packets corre- 
sponds to a situation where homogeneous line broadening is dominant, so that each neutrino 
wave packet is broadened because the energy of the emission line, Es, changes during the 
emission process. In contrast, for mostly inhomogeneous line broadening, each individual 
neutrino can be approximately described by a plane wave because it is emitted with an 
extremely small energy spread (which is ultimately determined by subdominant homoge- 
neous solid state effects, by the natural width, and by the Heisenberg principle). Different 
neutrinos, however, are emitted with different energies which depend, for example, on the 
proximity of the emitting atom to crystal impurities and lattice defects. 

To end this section, let us give a simpler and more useful form of eq. fl2^ . obtained by 
neglecting the localization and coherence terms and considering the two-flavor approxima- 
tion, with an effective mixing angle 6, an effective mass squared difference Am^, and an 
average absolute neutrino mass fh [l|: 
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3 Mossbauer neutrinos in quantum mechanics: Lorentzian wave 
packets 

Let us now discuss how oscillations of Mossbauer neutrinos can be understood in the 
framework of quantum mechanics. Since QM is unable to describe particle creation and 
destruction, we cannot directly include the production and detection processes into our 
formalism, as in QFT. Instead, we will first compute the probability for transitions be- 
tween the initial and final neutrino states, and then multiply this with the emitted flux and 
with the absorption cross section to obtain the overall event rate F. We will describe the 
propagating neutrino as a superposition of three wave packets, one for each mass eigen- 



state |32|, l33|, l3J, l35|, l36|, l37l|. As we have discussed above, such a description corresponds 
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to Mossbauer neutrinos in the regime of homogeneous hne broadening, while the case of 
inhomogeneous broadening would be more naturally implemented by considering an en- 



semble of many plane wave neutrinos in the density matrix approach [32|]. However, since 
homogeneous and inhomogeneous line broadening cannot be distinguished experimentally, 
it is sufficient to focus on one of the two cases. We use the wave packet picture because it 
provides insights into the evolution of each single neutrino, which we find useful to better 
understand the localization and coherence conditions that will emerge. 

Unlike most other authors, who use wave packets with a Gaussian shape, we will use wave 
packets with a Lorentzian momentum distribution because it is known from the classical 
Mossbauer effect that homogeneous and inhomogeneous line broadening mechanisms lead 



to a Lorentzian energy spread [38|, |39[ . The momentum space wave function for the electron 
antineutrino produced in decay thus has the form 



iPli^esit)) = ^Y1 "^"^ ,^ exp [ - zE,t\ 

Ns P-Pjs + i-ls/^ 



u,) . (32) 



The index 5* indicates that this state is produced in the neutrino source, and the normaliza- 

1 /2 

tion factor is Ns = ( J2j l^^ejP l/j5p) • Similarly, the detection process can be described 
as a projection of luesit)) onto a state Iucd) with the momentum space representation 



1 /2 

and the normalization factor = ( |f/ejp l/jSp) • In the above expressions, pjs, Pjd 
are the central momenta of the wave packets, Ej = [p^ + m|)^/^, and 75, 7/) are the wave 
packet widths. Moreover, we have introduced phenomenological fudge factors fjs-, fjo that 
will be motivated and discussed below. 

For Mossbauer neutrinos, 75 and 7£) are of the order of the energy uncertainty associated 
with the emission and detection processes, which is of order 10~^^ eV. The much larger 
momentum uncertainties of the source and the detector do not play a role because the 
neutrino is on-shell, so that by virtue of the relativistic energy-momentum relation the 
momentum uncertainty of the neutrino cannot be larger than its energy uncertainty. (Of 
course, the momenta associated with the different mass eigenstates have to differ by much 
more than 10~^^ eV in order to ensure energy-momentum conservation in the production 
and detection processes.) 



Note that lueo) is time-independent (on this point, we disagree with ref. |37|, where the 
detection operator \i'f3D){i'i3D\ is assumed to be not a time-independent but only a time- 
averaged quantity); on the other hand, a factor exp[—ipL] is required to center the wave 
packet around x = L. 

The phenomenological fudge factors fjs and fjc can be used to describe a possible mass 
dependence of the neutrino production and detection amplitudes. For example, we have 
seen in the previous section that the Lamb- Mossbauer factor depends on mj, so that the 
production and absorption of the lighter neutrino mass eigenstates is slightly suppressed 
compared to the production and absorption of the heavier ones. This can be viewed as a 
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slight dynamical reduction of neutrino mixing. Let us stress that fjs and fjn cannot be 
determined in the QM approach, and have to be put in by hand. We will choose 



fjs = exp 



fjD = exp 



(34) 



(with aps, <7pD, and E defined as in eqs. and (l20l) . respectively) in order to ultimately 
reproduce the correct Lamb-Mossbauer factor. 

The amplitude for the transition \iyes{t)) l^^eZ)) is given by 



A{t,L) = {UeDWesit)) = dp {VeD\v) {p\^es{t)) ■ 



(35) 



To be able to evaluate this integral, we make use of the smallness of 75 and 7d, and expand 
Ej around the average momentum pj = [pjs + Pjo)/'^, which gives 



E. 



'j = sjp^ + m] t Ejt + Vjt{p - Pj) 



with the definitions 



and 



Pj 



p] + 



(36) 



(37) 



This approximation corresponds to neglecting dispersion (wave packet spreading), which is a 



second-order effect j40|. Eq. (l36l) is a good approximation as long as {p — Pj)/Ej ^ E"^ / 



for all p within the peak regions of the source and detector wave packets. We can now 
compute v4.(t, L), and obtain 



NsNd 



X] \Uej\^ fjsfjD 



-Z75 



PjS - PjD -i{ls + Id)/^ 



exp [ — iEjt + iVjPjt\ 



■ jexp (^ipjs + y) ~ ^j^) ^(^i^ - L)+ exp (^ipjo - ^ j (L - vjt) 9{-Vjt + L) | , 

(38) 

where 9 denotes the Heaviside step function. 

The next step is to compute the transition probability for the process \i^es) li^eo), 
defined by 



1 

V{L) = - / dtA*{t,L)A{t,L). 

J- J-T/2 



(39) 



Here, the incoherent averaging over the running time T of the experiment refiects the fact 
that we do not precisely know at which point in time the production and detection reactions 
take place. (The detection time is, of course, implicitly constrained by the fact that the 
neutrino wave packet has sizeable overlap with the detector only during a very short time 
interval.) Physically, V{L) gives the time-averaged probability that a neutrino prepared 
in the state |i^eS'(0)) at t = is detected as \i'eD) at a later time. Note that 'P{L) is 
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not a z/g survival probability in the usual sense because in general V{L) |^„2 7^ 1. In 

jk 

particular, irrespective of the neutrino mixing parameters, V{L) can only be sizeable if the 
wave packets \i^es{t)) and Iued) have sufficient overlap in momentum space. This is precisely 
the Mossbauer resonance condition. 

The experimentally observable event rate F is obtained by multiplying V{L) with the 
Mossbauer neutrino emission rate F^^, the Mossbauer neutrino detection cross section cr^^, 
and the geometrical flux suppression factor 1 / AirL"^ : 

F = ^Fo^^P(L)a^^^ (40) 

- ^(rol5 ^ |f/e,f|/.5p) ViL) {B,Yr>^ Yl \Ue^V^n\') ■ (41) 
J 3 

The parenthesized expressions for F^^ and cx^^ have to be derived in the QFT formalism 
discussed in sec. |2]and ref. [ij. It is impossible to derive them in QM because they describe 
particle creation and annihilation processes. Note that we are here using the cross section 
for the limiting case of an inflnitely sharp Mossbauer resonance — hence the factor T/27r, 
which should be understood as an approximate 5-peak of the form 

5(0) ~ lim / rfte'(^-^^>°)* = — . (42) 



The effect of line broadening is already accounted for by the fact that A(t, L) is suppressed 
if \pjs - Pjd\ > (75 + 1d)/2 (cf. eq. dMD). 

Evaluation of F requires splitting the time integral in eq. fj39l) into three separate in- 
tegrals with integration domains {—00, L/vk], {L/vk, L/vj), [L/vj,oo) for mj > rrik, and 
{—00, L/vj], [L/vj, L/vk], [L/vk,oo) for nij < rrik- (It is justified to replace the integration 
boundaries ±T/2 from eq. ( l39l) by infinity here because the overlap of the wave packets 
l^esit)) and lueo) decreases exponentially at large T, when the neutrino has long passed the 
detector.) We will only show how to evaluate one of the above integrals, since the others 
are similar. Consider 



pL/vj 

Jjk = / dt exp 

J L/vk 



i{Ej - Ek)t + i{vjPj - VkPk)t - i{vjPjD - VkPks)t 

(43) 



+ ^ilDVj - -isVk)t + i{pjD - Pks)L -^{iD- ls)L 



for rrtj > rrik. We use the approximation of ultrarelativistic neutrinos {rrij ^ Ej), which 
suggests the expansions 

2 2 

TTl Tfl 

p,s ^ Es,o - (1 - e^)^ , PjD ^ En,o - (1 - ^d)77t;^ , (44) 

from which it follows that 

9 2 9 

J ^ ^ E ^ 2E^ ^ 
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where 



E 



5,0 



E 



D,0) 



and 



2 



1-e 



D 



E. 



5,0 



E 



D,0 



(46) 



In these expressions, Es^ and -E_d,o are the mean energies for the case of massless neutrinos, 
and C,s, C,D are constant parameters determined by the properties of the source and the 
detector, respectively. These parameters can be calculated only in an explicit treatment of 
the neutrino production and detection processes. For conventional neutrino sources, C,s and 

are of 0{1), but for Mossbauer neutrinos, the energies associated with different neutrino 
mass eigenstates have to coincide within the line widths 75 and •jd, so that C,s, (,d and ^ 
must be extremely small in this case. Indeed, we will see below that, in order to reproduce 
our QFT result (El), we have to take = = 0. 

Plugging (145|) into (l43l) . neglecting terms containing the small product Am|^(i?5^o — 
Eo,o)/E'^ and, in the denominator, also neglecting terms of order 'ym'j/E'^, we obtain 



J, 



jk 



AS) 



AD) 

^jk 



|(7D - 75) + ^(^5,0 - Ed,o) - t^Am%/2E 
with the oscillation and coherence terms abbreviated as 



(47) 



AB) 
^jk 



exp 



— z 



Am%L 



= exp 



- 27rz- 



2E 
L 



\^m%\lBL 
4^2 



r osc 
^jk 



T coh 

^BJk. 



for B = {S,D}. Note that the A^-^^ are identical to the quantities of the same name 
defined in eq. ( 122]) . up to the replacement of Eb,o by E, which leads to corrections of 



0{Am'^i^{Es^o — Eofl)/ E"^). Since we have neglected terms of this order in the derivation of 
we should for consistency also neglect them here. The full expression for F is 



47rL2 



1 

'2^ 



l^lf/efel^ exp 



2^2 



m 



2l 



ISlD 



E. 



5,0 



E 



E. 



5,0 



E 



Dfl 



m,; 



ml 



A 



is) 

jk 



. 75 + i^i 



' 2Es,o 



is 



2E, 



5,0 

AD) 
^jk 



2Ed,o 



2E 



D,0 



Id - i^D 



+ 



i{ls + Id) 



A 



(D) 
jk 



UlD - 75) ± i{Esfi - Ed,o) - 



2E 



(49) 



In the last term, the upper sign applies to the case Aim}-j^ > 0, while the lower one is for 
Am^^ < 0. As discussed above, ^5 and are very small for Mossbauer neutrinos. If we 
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neglect them completely, T simplifies to 



47rL2 



2tt 



ej I 



Is +lDf . (5) 



jk ) 



\\Jekf exp 



1 (^i-? - 4"') [(^5,0 - i?Ao)(7s - 7d) ± 



(E5,o-^D,o)' + i(75 + 7D)^ 



; (75+7p)^ 



S,0 



• (50) 



This equation is identical to our QFT result, eq. (IMl) (within the approximations made in 
the two approaches). In particular, we find the same oscillation, coherence, and resonance 
terms. 



4 Discussion and conclusions 

Let us now summarize and discuss our results. In the first part of this paper, we have 
used quantum field theoretical techniques to derive the rate F of Mossbauer neutrino emis- 
sion, propagation, and absorption (eq. (12^ ). For the first time, we have explicitly included 
the effect of homogeneous line broadening due to fluctuating electromagnetic flelds in the 
solid state crystals forming the source and the detector. We have conflrmed the expecta- 
tion from ref. [l| that the resulting formula for F agrees precisely with the one obtained 
in [l[ for the case of inhomogeneous line broadening caused by crystal defects and impuri- 
ties. In particular, we have conflrmed that, also for homogeneous line broadening, F has a 
Breit-Wigner-like resonance structure, and contains oscillation, localization, and coherence 
exponentials. Moreover, our formula accounts for the suppression of recoilless emission and 
absorption processes compared to their non-recoilless counterparts through a generalized 
Lamb-Mossbauer factor. We have also noted that in realistic experiments the localization 
and decoherence terms are irrelevant and may be set equal to unity. The localization term 
enforces the condition that the quantum mechanical delocalization of the neutrino source 
and detector have to be small compared to the oscillation lengths for oscillations to take 
place, a condition that is easily fulfllled in any oscillation experiment. The decoherence 
term, on the other hand, accounts for the possibility of wave packet separation due to the 
different group velocities associated with different neutrino mass eigenstates, but also this 
does not happen in terrestrial experiments. 

We have then proceeded to a derivation of F in a quantum mechanical approach, in 
which the neutrino is described by a Lorentzian wave packet of the form (1321) . We have 
arrived at eq. (l5Up . which coincides with the QFT result (1241) . However, since the neutrino 
production and detection processes, which involve particle creation and annihilation, cannot 
be described in QM, the Mossbauer neutrino production rate as well as the detection cross 
section had to be put in by hand. Also, the properties of the neutrino wave packets (shape, 
width, central momenta) had to be chosen in an ad hoc way instead of emerging naturally 
from the formalism or being related to properties of the source and the detector. Once 
the appropriate choices for these parameters are made, the Breit-Wigner-shaped resonance 
factor as well as the oscillation and decoherence terms can be derived. The correct Lamb- 
Mossbauer and localization factors are obtained only if suitably chosen phenomenological 
weighting factors /j^, /j^) for the different neutrino mass eigenstates are introduced in the 
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neutrino wave function to account for the tiny dependence of the emission and absorption 
probabihties on the neutrino mass. 

As expected, F factorizes into the emitted neutrino flux, a transition probabihty V{L), 
and the detection cross section. While in the QM approach, this property is introduced 
as an assumption in eq. fHOl) . it emerges naturally in QFT. The reason is that for large 
propagation distance L off-shell effects become negligible, and according to the Grimus- 
Stockinger theorem eq. ^ the propagator then reduces to the exponential phase factor 
exp[ipL) (with p being the modulus of the neutrino momentum), which is also used in QM 
to describe the spatial evolution of particles. 

In conclusion, we have shown that the QM approach to Mossbauer neutrino oscillations, 
in which the production, propagation, and detection processes are treated separately, is 
able to reproduce the results obtained in the QFT approach, in which these processes are a 
priori considered as a single entity and their factorization emerges as a result. In general, the 
framework of QFT is significantly more robust because it does not require any assumptions 
on the neutrino wave function, whose parameters are instead automatically determined 
from the much less ambiguous properties of the neutrino source and the detector. For 
example, homogeneous and inhomogeneous line broadening are easy to implement in QFT 
(see sec. |2]and ref. while in QM, they have to be accounted for by choosing appropriate 
wave packet widths. Also, the emission rate, the detection cross section, and the Lamb- 
Mossbauer factor cannot be predicted in QM and have to be put in by hand. On the other 
hand, the QM approach can give a better physical understanding of the origin of oscillation, 
decoherence, and resonance phenomena once all free parameters are chosen appropriately, 
e.g. by matching with the QFT result. 
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